Failure of the displaced-squeezed state for spin-boson models in the thermodynamic 

limit 
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We present an analysis of a variational coherent-squeezed state that has been discussed in the 
literature as a potential ground state for the spin-boson model. We show that when the system-size 
scaling of the spin-bath coupling is included properly, all squeezing effects and non-universal physics 
vanish in the thermodynamic limit. We also present finite-size corrections to the renormalisation of 
the spin's coherence, showing that squeezing effects are also absent to leading order in the inverse 
bath-size. 



I. INTRODUCTION 

The spin-boson model is one of the most important 
theoretical models for studying dissipation and decoher- 
ence in quantum systems, and has been applied to numer- 
ous systems in chemistry biology 0, and the emerging 
fields of quantum computation and quantum devices 3 . 
The model is simple, comprising a coherent two level sys- 
tem (TLS) that is coupled to a large bath of harmonic 
oscillators, yet despite this simplicity, the model is excep- 
tionally rich and continues to provide new insights into 
open many-body quantum physics. 

One particularly important issue, especially for practi- 
cal quantum devices, is the robustness of the coherence of 
the TLS in equilibrium with the bath. Numerous studies 
on this have revealed some extremely interesting physics, 
including quantum phase transitions between coherent 
and incoherent TLS states, and unusual tem per ature de- 
pendence of the cohrencejl S !, 0, H S 01, El ■ 

Amongst these studies, the use of trial wave functions 
has provided some very intuitive insights into the prob- 
lem. There are two commonly discussed states which are 
both based on, and improve on, the adiabatic approxi- 
mation discussed by Leggett et al ll| . The first was that 



of Silbey and Harris, who modified the adiabatic approx- 
imation to allow for simple non-adiabatic responses from 
the low frequency modes in the problem^. This state 
was shown to correctly describe the coherent-incoherent 
transition for Ohmic coupling^, and has recently been 
used to study the transition in the sub-Ohmic system, 
giving results in excellent agreement with those obtained 
by other methods @, 0, @, Gj| ■ 

The second state is known as the displaced-squeezed 
state (DSS) and was proposed by Chen, Zhang, and 
Wong[12|. This state is based on two effects; adiabatic 
displacement of the bath modes, and spatial deformation 
of the oscillator wave functions. At zero temperature it 
has been shown that the energy of the DSS can be signif- 
icantly lower than the Silbey-Harris state for the case of 
strong coupling between one oscillator and the TLS[l~3l|. 



and this has also been claimed for strong coupling to a 
super-Ohmic continuum of bath modes[14|. 

This DSS also predicts coherent-incoherent transitions 
in Ohmic and sub-Ohmic systems [lH [l6j. but with sig- 
nificantly different critical properties than those found 
in the rest of the literature. Moreover, these critical 
properties depend on the explicit form of the TLS-Bath 
coupling, which is unusual as it has been shown us- 
ing path integral techniques that the spin-boson physics 
should be controlled solely by the spectral function of 
the bath [n|, G3- As the DSS* is claimed to be the most 
stable ground state for strong coupling, there has been 
much discussion of this breakdown of 'universality' as a 
function of coupling strength [Til, [f5j . a breakdown which 
we shall show does not exist thermodynamically. 

In this paper we use an effective Hamiltonian theory 
to show that when the system-size scaling of the TLS- 
bath coupling is included properly, all squeezing and non- 
universal features of the DSS vanish in the thermody- 
namic limit. The physics which emerges is simply that 
of the basic adiabatic approximation. Finite-size effects 
are then discussed, with no evidence for squeezing found 
to leading order in the inverse system-size. 



II. THE SPIN-BOSON MODEL 

We shall use a standard picture of the spin-boson sys- 
tem in which we consider the TLS as a particle tunnelling 
between the minima of a double well potential as shown 
in figHJ The Hamiltonian for the spin-boson model is 
given by[nj], 



-Ka, 
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9l(ai + a}) + 



(1) 



where a XtZ are Pauli matrices, ai,aj are the bosonic cre- 
ation and annihilation operators for the bath modes, gi 
is the coupling between the TLS and the Ith bath mode, 
lui is the frequency of mode I, and K is the tunnelling 
matrix element between localised states. We have cho- 
sen the TLS basis such that the state localised in the left 
hand well is | f ) and | J.) is the right well state. 
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FIG. 1: The Coherent Two- level system is modeled as particle 
tunnelling in double well potential. Bath oscillators become 
correlated with the position of the particle, renormalising the 
matrix element for tunnelling due to the poor overlap of their 
displaced wave functions. 



Following Leggett et al we assume that the bath is 
truly macroscopic and that the modes can be treated in 
the continuum limit with a smooth density of states up 
to some cut-off frequency u> c [H| ■ The physics of the spin- 
boson model is then normally determined solely by the 
spectral function of the bath J(u>), 



aw s s lo s 9(uj c - lj). 



(2) 
(3) 



The RHS of eq.Q is a phenomenological power-law 
where a is a dimensionless coupling strength, uj s is a 
typical frequency scale of the bath, and s is the exponent 
of the frequency dependence. The ultra-violet cut-off lo c 
for the spectral function will be taken to be the largest 
energy scale in the problem, much larger than K or u> s . 
This form for J(u>) can derived for many specific micro- 
scopic interactions, e.g. phonons, E.M modes, and many 
other physical examples can be found in the literature [J]. 

The dynamical and thermodynamical behaviour of the 
TLS depends critically on the coupling strength and the 
exponent s. As a result, three types of bath arc distin- 
guished in the literature according to their value of s, the 
super-Ohmic bath (s > 1), Ohmic bath (s = 1), and the 
sub-Ohmic bath (s < 



These ground states can be understood intuitively; the lo- 
calised TLS simply creates a static force which displaces 
all the bath modes, the displacement being described by 
the action of the shift operator Ua on the vacuum of all 
bath modes |0). When K is finite the TLS can tunnel 
between the wells, and no exact solution for the problem 
is known. However, bath modes with frequencies much 
higher than the tunnelling rate 2K can respond almost 
instantaneously to the slow tunnelling motion and can be 
eliminated using an adiabatic approximation 

In the zeroth-order adiabatic approximation, where all 
modes are assumed to follow the TLS instantaneously, 
the tunnelling leads to a coherent ground state \g.s) = 
y=(|+) + |— )), characterised by a renormalised coherent 

level splitting K = (+\K<j x \—) between the tunnelling 
states |+) and |— ). Physically, the tunnelling probability 
is reduced by the TLS-bath correlation, as the overlap 
between the states |+) and |— ) is suppressed by the rel- 
ative displacement of the oscillators as shown in figffl 
Calculating K explicitly, and using the definition of the 
spectral function to write the sum over bath modes as an 
integral, we find that, 



K exp 



(6) 



From Eq. ([6]) and Eq. ([3]) we can see that K is finite for 
super-Ohmic baths, and zero for Ohmic and sub-Ohmic 
baths due to the infra-red divergence of the integral in 
Eq. ([HI). Thus the zero-order adiabatic approximation 
predicts that the TLS is always localised in an Ohmic 
or sub-Ohmic environment at T = K, a result which 
is known to be incorrect @, 0, H, @- This failure can 
be traced to the mistreatment of the slow modes in the 
problem, for which the adiabatic approximation is clearly 
not valid. 

The DSS of Chen, Zhang and Wu[l2| attempts to im- 
prove on the adiabatic approximation by also allowing for 
spatial distortion of the oscillator wave functions as they 
adiabatically follow the TLS. As is discussed in Ref. [121 ]. 
the distortion of the oscillator wave functions can be de- 
scribed using the generators of bosonic squeezed states, 
minimum uncertainty states commonly used in quantum 
optics The effective displaced-squeezed ground state 



they propose is given by \ips) 
U s is given by, 



U A U S 



(IT)+li))|Q) 

V2 



where 



III. ADIABATIC APPROXIMATION AND THE 
DISPLACED-SQUEEZED STATE 

In the absence of tunnelling the TLS will be localised in 
one of the wells, and the spin-boson Hamiltonian can be 
solved exactly. The two degenerate ground states |+), |— ) 
are given by, 

|+) = u A \W), \-) = u A \i)\o), (4) 

U A = cxp[-(7 z ^giuf 1 (a i - a])}. (5) 
i 



l/ s = exp[-£ 7( (a?-a| 2 )]. (7) 
/ 

The parameters 7/ describe the amount of spatial distor- 
tion of the oscillator wavefunction, or alternatively, the 
squeezing of the bosonic quadrature operators [l8j|. Once 
determined, the squeezing parameters modify the renor- 
malisation of the tunnelling matrix clement by altering 
the overlap integral of the adiabatically displaced oscil- 
lator states. This modification can lead to qualitatively 
different physical results, including the possibility of fi- 
nite K for Ohmic and sub-Ohmic systems. 
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IV. VARIATIONAL METHOD AND THE 
EFFECTIVE HAMILTONIAN 

We now use an adapted version of the variational 
method of Silbey and Harris to calculate the effective 
tunnelling matrix clement of the TLS in the displaccd- 
squeezed state at T — K[!5]. First we make a canonical 
transformation to generate the Hamiltonian in the basis 
of the displaced-squeezed state, H = UAUsH s i,Ug 1 U^ 1 . 

H = —K + <7 + — K-O- + u)\ sinh(27;) 

I 

- ^ff^r 1 +^^aJa ; (cosh 2 (2 7; ) +sinh 2 (2 7; )) 



+ ^tuiiaf + a\ 2 ) cosh(27;) sinh(2 7 z), 



where the operators K+, K_ obey, 



K 4 



K* 



K exp 



(8) 
(9) 

(10) 



In the approximations discussed previously, the inter- 
acting bath-TLS system can be described by an effective 
non-interacting Hamiltonian characterised by the renor- 
malised tunnelling matrix element K. We derive this 
form of Hamiltonian, as a mean-field approximation to 
H , by introducing the expectation value of the tunnelling 
operators K = (0|.K+|0) = (0|^_|0), which has the ex- 
plicit form, 



K = K exp 



-2$>V 2 e-^ 



(11) 



Adding and subtracting Ka x to H , we then write the 
Hamiltonian as H = Hq + V, 



Ha 



~Ka x + ^2 ui sinh(27i) - ^ gfu 
i i 

y^^q]'a;[cosh 2 (27;) + sinh 2 (27f)] 



(12) 

V = H-H Q . (13) 

Following Silbey and Harris [||, we now compute the 
Bogoliubov-Fe ynm an upper bound on the free energy of 
the system Aq [lj| . Bogoliubov's theorem states that the 
true free energy A of the model is related to Ab bv[2fj|. 



A < A, 



(14) 



At 



-/T 1 lnTrexp(-/3 J ff ) + <V>fl- . (15) 



The angular brackets denote the thermal expectation 
value calculated with respect to Hq, and computing the 
trace using the simple eigenstates of H , we find that 
(V)h = 0. The free energy bound at T — K is thus, 



-K + sinh ( 2 7i) - J2 Vi 



V 1 - 



(16) 



We note here that, like all previous studies in the lit- 
erature, this variational method does not consider the 
possibility that the separation of the Hamiltonian into 
Hq + V may lead to divergent fluctuations if higher or- 
der perturbations are calculated i.e. that the variational 
state is an inappropriate starting point for analysis of this 
problem. This matter will be discussed in a forthcoming 
study of these variational methods. 

The variational parameters 7; are determined by min- 
imising the energy of the system w.r.t the set {72}. Wc 
find that 7; satisfies, 



8Kgf 



(17) 



and substituting this into Eq. pTj) . wc find that K obeys 
the self-consistent equation, 



K = K exp 




(18) 



This self-consistent equation must be solved to find 
the effective tunnelling matrix element of the mean field 
effective Hamiltonian Hq. For a finite K the ground state 
is a delocalised, coherent superposition state, whilst for 
K = the TLS is incoherent and localised in one of 
classical well states. Due to the self-consistent nature of 
Eq. p^|) . the function K may vanish at a discrete value 
of the coupling strength a, signifying a type of phase 
transition between coherent and incoherent phases 0, 0, 
110. 

The presence of gf in two places in Eq. (TT5)) means that 
the renormalisation of K is determined by both the spec- 
tral function and gf. In the literature [H, [HI, EH, <?z is 

taken to have a general power-law form gi = 9a (^^j , 
where go is a constant. However, this form for the cou- 
pling constants cannot be valid as it doesn't not take 
into account the scaling of these constants with system 
size. From Eq. (f5|) we see that J (to) is the product of the 
square of the coupling constant at lji = lu and the density 
of states per unit frequency. As the density of states per 
unit frequency is proportional to A/", where Af is the total 
number of bath oscillators, the microscopic coupling con- 
stants have to scale as N~ 2 to ensure that the spectral 
function is well defined in the thermodynamic limit. This 
point is made very clear in the origin review of Leggett 
et al[H|. 

The non-scaling form for gi used in the DSS literature 
may be related to the continuous coupling function g{oS) 
that can be defined from the spectral function and bath 
density of states, as discussed by WeissQ- The spectral 
function can be written as J(u>) = 5 2 (w)2?(w)[3|, where 
T>(u>) is the bath density of states per unit volume. This 
coupling function g(to) can take the non-scaling form used 
in the DSS literature, as the intrinsic scaling of the micro- 
scopic couplings gi has been absorbed into the definition 
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of the density of states per unit volume. However, the 
quantities that appear in the self-consistent equation Eq. 
(fl~8|) are the actual microscopic coupling constants from 
the spin boson Hamiltonian, and these must include the 
system-size scaling. As the dependence of AT on Af is es- 
sential to our discussion, we assume a general power-law 
form for the couplings given by, 



9i 



(19) 



where n an exponent greater than zeroQ. 

We now use the definition of the spectral function given 
in Eq.© to write the sum in Eq. fTg)) as an integral over 
a continuous distribution of bath modes. In this limit, 
the self-consistent equation for K takes the form, 



A" 



K exp 



-2 



J(lo)u> 2 dio 



1 + %auyt n Kw n ^M- 1 



K = A" exp -I(K) 



(20) 



We then write the integral I{K) in dimcnsionlcss form, 

w s ~ 2 dio 



I{K) = ato 



I(K) = a 



l-s 



(l + SaKco^co^Af- 1 

!-s pi „s-2 



X s dx 



1 + A(K)x 



n — 3 



i * 



(21) 



where we have defined the dimcnsionlcss numbers x = — 

and A(K) = Comparing Eq.([2T]) with Eq.© we 

can see that the inclusion of squeezing effects introduces 
a square root factor into the integrand which can have a 
profound effect on the renormalisation of K. The squeez- 
ing provides an effective infra-red cut-off to the integral, 
and thus prevents the infra-red divergence if n < 2s + 1 
and A(AT) is finite. Squeezing effects can therefore po- 
tentially lead to coherent ground states for Ohmic and 
sub-Ohmic baths, as well as sharp coherent-incoherent 
transitions, as is described for the Ohmic case in refs. 

However as there can only be a true infra-red diver- 
gence in the thermodynamic limit, we must consider 
what happens to the solutions of Eq. ([20|) as Af — > oo. 
In this limit we see that A — > 0, and thus the inte- 
gral in Eq. (|21[) reverts to the form given by the adia- 
batic approximation [H [ll|. Therefore for Ohmic and 
sub-Ohmic baths, the DSS always leads to an incoher- 
ent thermodynamic ground state, which from Eq. (|17p 
leads self-consistcntly to the vanishing of the squeezing 
parameters ji. Interestingly, a recent variational study 



of squeezing effects in a weakly interacting Bose-Einstein 
condensate also found that the squeezing vanishes in the 
thermodynamic limit. [22l|. 

We can also show that this conclusion is independent 
of when we choose to take the thermodynamic limit by 
solving the self-consistent equation at finite Af and then 
sending Af — > oo. At finite Af there can be no phase 
transition as there will always be an infra-red cut-off lojr 
arising from the finite-size of the system. This frequency 
natural scales inversely with the linear dimensions of the 
system, and can be written loir = J\f~3ui 0y where u> 
is some system dependent frequency which depends on 
the density of oscillators Af/V. The squeezing effects 

provide an effective infra-red cut-off at a frequency of 

i i 

approximately w c A 3 -™ oc Af 3 ~" , which vanishes faster 
than loir as Af becomes very large. Thus in the limit 
of large Af, the squeezing factor can be ignored in the 
integrand of Eq. (|21|) . and the renomalisation is given 
by the adiabatic formula of Eq.©, but with an infra-red 
cut-off at loir. The renormalised matrix element is then 
given by, 



K = K exp 



sAfi 



1 - s 



UJQ 



(22) 



for the case of sub-Ohmic coupling, and 



K = A 



too 



J c Af3 



(23) 



for Ohmic coupling. 

As expected, for finite Af the TLS retains it's coher- 
ence for arbitrary coupling strength, with the coherence 
vanishing in the thermodynamic limit. Due to the fact 
that the infra-red divergence is prevented by lojr rather 
than lo c As^ , the coherence does not depend on the TLS- 
bath coupling, although these finite-size corrections are 
non-universal as they depend on the bath frequency too. 

Although we find no evidence for squeezing even at 
finite Af, squeezing effects can be very important for 
coupling to a single mode[l3[, or possibly a few discrete 
modes. For instance, such squeezing effects have recently 
been predicted for a TLS coupled to a single nanomc- 
chanical oscillator [23l|. 

In conclusion, when one uses the correct system-size 
scaling of the TLS-bath couplings, analysis of the DSS 
shows that all the squeezing and non-universal effects 
vanish in the thermodynamical limit. In thermodynamic 
equilibrium, this variational state thus reproduces the 
results of the basic adiabatic approximation, and fails to 
capture the interesting transition physics of Ohmic and 
sub-Ohmic systems. Squeezing effects arc also found to 
be absent when one considers the leading order finite-size 
corrections around the zero-order adiabatic state. 
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